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Abstract 

Simulations  of  the  performance  of  multi-spectral  sensors  or  multiple 
sensors  with  data  fusion  require  phenomenological  consistency  in  the  target 
and  background  signatures  presented  to  the  sensors.  This  paper  describes 
the  incorporation  of  ultraviolet  signatures  of  targets  and  backgrounds  into 
an  engagement  simulation  framework.  This  paper  emphasizes  ultraviolet 
phenomenology  and  the  signatures  of  aircraft  compared  to  fires,  lights,  in¬ 
dustrial  sources  and  other  high  energy  sources.  In  particular,  variation  of  the 
sources  as  a  function  of  time,  wavelength  and  aspect  angle  is  discussed,  as 
these  variations  may  be  critical  to  the  discrimination  of  targets  from  natural 
and  man-made  background  sources.  The  organization  of  multiple  data  sets 
representing  the  same  object  in  different  spectral  regions  is  also  described 
for  use  in  simulations. 
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1  Introduction 


Traditionally  GTSIMS  missile  flyout  simulations  have  consisted  of  a  mis¬ 
sile,  a  target  and  a  background.  Backgrounds  were  composed  primarily  of  natural 
features  such  as  sky,  water,  rocks,  roads,  fields  and  trees.  These  features  are  gen¬ 
erally  adequate  since  most  infrared  missiles  have  small  field  of  views,  short  flight 
times  and  upward  launch  geometries,  making  the  possibility  of  other  man  made 
objects  in  the  background  low. 

Under  the  ATIRCM/CMWS  program  we  have  been  refining  the  capability 
to  test  optical  missile  warning  receivers  (MWRs)  as  players  in  closed  loop  sim¬ 
ulations.  Optical  missile  warning  receivers  employ  sensors  with  a  up  to  a  An  sr 
field  of  regard,  and  these  devices  scan  this  field  of  regard  for  sources  that  match 
the  characteristics  of  threat  missiles.  The  large  fields  of  regard  of  these  systems 
along  with  their  need  to  watch  the  ground  for  long  periods  of  time  translate  to  a 
high  probability  that  man-made  objects  other  than  the  threat  will  be  visible  to  the 
MWR  during  a  simulation. 

Threat  missiles  typically  have  energetic  plumes  that  are  brighter  than  most 
natural  background  objects.  However,  a  number  of  common  objects  such  as  lights, 
fires,  aircraft  plumes,  welders,  pyrotechnic  flares,  bombs  and  other  weapons  can  a 
be  least  as  energetic  as  missile  plumes.  As  part  of  the  background,  these  sources 
can  have  two  major  effects  on  the  systems  being  tested.  First,  they  can  cause  false 
positives  which  reduce  pilot  confidence  in  the  system.  Second,  as  clutter  objects, 
they  can  conceal  the  true  threat. 

In  order  to  properly  test  systems  in  the  presence  of  these  sources,  a  database 
of  sources  was  constructed  from  measured  data  and  predictive  codes.  Framework 
was  added  to  the  simulation  so  that  many  of  these  sources  can  be  easily  added  to 
the  simulation. 

This  paper  discusses  flexibility  provided  by  the  model  framework,  lists  the 
background  and  target  sources  available  in  the  simulation  and  discusses  efforts  to 
model  specific  sources  in  the  ultra-violet  (UV)  region  of  the  spectrum. 


2  Object  Descriptions 

Any  object  in  the  simulation  that  emits  or  interacts  with  radiation  from  UV 
through  IR  wavelengths  is  summarized  by  a  visible  object  file  which  tells  the  sim¬ 
ulation  how  to  present  the  object  to  a  different  sensors  that  are  in  the  simulation. 
The  different  classes  of  sensors  are  mainly  defined  by  their  wavelength  character- 


istics  such  as  UV,  near  IR,  mid  IR  etc.  The  classes  also  separate  sensors  that  are 
likely  to  resolve  the  target  with  many  pixels  and  process  the  resulting  spatial  pat¬ 
tern  from  sensors  that  are  never  likely  to  resolve  the  target.  Examples  of  resolved 
and  unresolved  sources  in  the  simulation  are  shown  in  figure  1 . 


Figure  1:  An  F-16  in  full  afterburner  shown  as  a  resolved  infrared  source  and  an 
unresolved  UV  source.  Range  to  target  is  approximately  IKm.  IR  field  of  view 
is  1.5  degrees  and  UV  field  of  view  is  100  degrees.  LOWTRAN  7  was  used  for 
IR  atmospheric  effects  and  OSIC  7.0  was  used  for  UV  effects.  Images  have  been 
processed  to  show  relevant  features. 

Resolved  objects  require  inputs  to  the  visible  and  infrared  scene  Tenderer, 
GTVISIT.  These  inputs  include  faceted  target  geometry  files,  surface  tempera¬ 
tures,  surface  optical  properties,  plume  flowfields,  and  matching  atmospheric  pa¬ 
rameters.  The  GTVISIT  software  is  described  in  [1]. 

For  unresolved  objects  one  of  the  most  important  cues  is  the  signature 
variation  verses  time.  To  a  sensor,  an  object’s  temporal  variation  depends  on  its 
inherent  temporal  signature,  its  change  in  location  with  respect  to  the  sensor,  its 
change  in  orientation,  and  changes  in  spectral  output  that  change  the  atmospheric 
transmittance  from  the  source  to  the  sensor.  For  short  wavelengths,  the  scatter 
pattern  that  reaches  the  sensor  can  also  be  an  important  detail  and  it  is  dependent 
on  the  object’s  spatial  distribution  as  well  as  atmospheric  properties. 

In  order  to  represent  each  of  these  aspects,  the  radiant  intensity  of  the  point 
object,  /  is  dependent  on  two  aspect  angles,  0  and  4>,  time,  t,  altitude,  a,  and 
velocity  v.  In  the  simulation  the  radiant  intensity  of  a  source  can  be  expressed  as 
the  product  of  multipliers  shown  below 


(1) 


7(A,  6,  0,  t,  St,  a,  v)  =  Sf()i,  8)i)J(t)A(6,  cp) P {a)V (v) 

where  A{0,(p)  is  an  aspect  angle  function  that  defines  the  relative  magnitude  of  the 
source  viewed  from  different  angles.  P{a)  is  a  multiplier  dependent  on  altitude 
and  V  (v)  is  a  multiplier  dependent  on  velocity.  The  spectral  multiplier,  S/(A,  SX) 
defines  the  portion  of  7(0  that  is  distributed  between  A.  and  A.  +  5A.. 


SfiX,SX) 


S(X)dX 
rXmax  S(X)dX 

^  Amin 


(2) 


where  5(A)  is  a  spectral  distribution  function  and  X^in  and  X^ax  are  the  minimum 
and  maximum  wavelengths  in  the  point  signature  definition. 

In  the  cases  where  the  signature  can  be  defined  as  a  sum  of  products, 
the  functions  in  equation  1  are  input  as  a  series  of  ASCII  tables.  The  ASCII 
tables  are  generally  small  and  loop  constructs  are  provided  to  play  back  peri¬ 
odic  sources.  When  there  is  coupling  between  the  functions,  and  entire  table  of 
7(A,  9,  0,  t,  St,  a,  v)  values  can  be  entered  directly  as  a  binary  table  which  is  en¬ 
coded  for  portability  across  machines. 


Class 

Specific  Examples 

Flares 

M206,  MJU7A/B,  MJUIO,  MJU8B,  MK46, 

MJU47,  MJU48,  XM212,  LUU2B 

Fires 

Diesel,  gasoline,  wood,  JP-8 

Lights 

High  intensity  discharge  and  quartz  tungsten  halogen 

Welders 

Stick  welders,  gas  tungsten  arc  welders 

Aircraft 

F-16  Pratt  and  Whitney  FI 00  and  GE  110  Engines,  FI 5 

Munitions 

m57c02,  m57c08,  m57cl0,  m831al,  m865e3, 
m934  mk82  mk83  mk84 

Table  1 :  List  of  point  sources  available  to  the  simulation. 


Sources  that  are  included  in  the  simulation  are  listed  in  table  1.  The  muni¬ 
tions  signatures  are  contributed  by  SciTec  and  are  visible  only  in  the  UV  portion 
of  the  simulation. 


3  Trajectories  and  Simulation  Geometry 

Another  cue  for  the  discrimination  of  unresolved  sources  is  motion.  Ob¬ 
jects  can  be  eliminated  from  consideration  if  predictions  of  their  velocities  are  too 
high  or  too  low.  Other  motion  cues  can  be  obtained  from  predictions  of  object 
direction  relative  to  the  sensor. 

In  the  simulation  object,  trajectories  originate  from  one  of  two  places. 
Missiles,  some  aircraft,  and  counter  measure  flare  trajectories  are  products  of  the 
closed  loop  simulations.  In  these  cases,  an  initialization  phase  communicates  in¬ 
formation  in  the  visible  object  file  to  executables  that  need  it  to  generate  imagery 
for  a  particular  sensor.  Once  the  initialization  is  complete,  position  and  orien¬ 
tation  updates  are  sent  as  they  are  required.  Other  background  objects  also  also 
associated  with  motion.  Examples  are  rotating  navigation  beacons,  illumination 
flares  like  the  LUU2B,  and  sources  on  missiles  or  other  aircraft  not  involved  in  the 
closed  loop  portion  of  the  simulation.  For  these  sources  a  simulation  object  was 
added  to  read  a  list  of  point  sources  and  broadcast  the  initialization  and  kinematic 
update  messages  for  each  source.  The  point  source  manager  uses  a  cubic  spline 
fit  to  trajectory  data  for  position  and  orientation  updates. 


Figure  2:  Billboard  with  Metal  Halide  Lamp  in  the  visible,  UV  and  UV  temporal 
trace.  One  damaged  lamp  and  its  reflection  is  visible  in  the  UV  image  collected 
by  the  MWIP  program  [3]. 

As  mentioned  above,  change  in  relative  positions  of  the  sensor  and  object 
changes  the  apparent  temporal  trace  seen  by  the  sensor.  This  modulation  primarily 
results  from  range  squared  fall  off  and  change  in  atmospheric  transmission  with 
range.  Another  source  of  trajectory  related  modulation  is  local  obscuration  of  the 
source.  Figure  2  demonstrates  a  common  example  of  this  phenomena.  A  damaged 


metal  halide  lamp  illuminates  a  billboard.  From  one  side  of  the  billboard,  the  lamp 
and  its  reflection  are  clearly  visibe,  but  neither  are  visible  from  the  opsite  side. 

The  GTVISIT  software  keeps  track  of  some  obstructions  between  sensors 
and  the  target,  but  this  technique  is  not  generally  applicable  for  missile  warning 
receiver  applications  because  the  large  area  visible  in  the  warning  receiver  field  of 
view  and  a  small  facet  size  needed  to  reduce  surface  error  near  the  source  require 
too  many  facets  for  practical  rendering  times.  For  geometries  involving  walls  or 
rectangular  enclosures,  a  utility  is  provided  to  compute  obstruction  and  diffuse 
reflection.  The  output  of  the  utility  is  an  aspect  angle  dependence  function  which 
can  be  added  to  the  point  signature  file.  Figure  3  demonstrates  a  fly  over  of  the 
from  an  the  back  of  a  simulated  billboard  illuminated  by  a  single  lamp.  When  the 
sensor  approaches  from  the  back,  it  never  directly  sees  the  lamp  or  the  reflection. 
As  a  result,  the  recorded  signal  decreases  as  the  sensor  apporaches  because  the 
billboard  blocks  a  larger  portion  of  the  scattered  light  from  the  source. 
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Figure  3:  Sensor  output  when  approaching  the  billboard  from  the  front  and  the 
back.  The  sensor  is  pointed  forward  and  passes  40  meters  above  the  source.  Initial 
range  is  500  meters. 


The  temporal  characteristics  of  the  lamp  used  in  this  example  come  from  ex¬ 
amination  of  Metal  Halide  and  Mercury  Vapor  lamps  measured  under  MWIPS 
and  by  the  CMWS  program.  Based  on  these  measurements,  a  fit  for  the  radiant 
intensity  of  the  lamp  is  defined 


J(t)  =  Jo 


sin 


*  8 


(3) 


Where  Jo  is  the  peak  intensity  of  the  lamp  output,  N  is  the  sample  rate 
in  samples/sec.  cod  =  (o/D  where  (o  is  the  power  line  frequency  in  Hz  and  D 


is  a  duty  cycle  factor.  3()  is  the  Dirac  delta  function,  and  i  <  N / 2.  The  tem¬ 
poral  distance  between  lamp  pulses  is  l/2co  and  the  width  of  the  pulses  at  zero 
amplitude  is  l/2cod.  In  some  cases,  the  resulting  time  sequence  is  filtered  with 
a  a  single  time  constant  filter  3.  The  peak  radiant  intensity  is  a  function  of  lamp 
wattage  and  construction.  Magnitudes  for  intensity  and  spectral  output  are  taken 
from  measurements  and  references  on  UV  light  output  such  as  [4]. 

Since  the  lamps  are  small,  they  are  also  treated  as  point  sources  in  the  in¬ 
frared.  The  primary  source  of  output  from  the  IR  lamps  was  assumed  to  be  emis¬ 
sion  from  the  glass  bulbs,  heated  fixtures  and  filaments  of  quartz  lamps.  Steady 
state  thermal  analysis  was  performed  to  predict  average  temperatures  of  lamps.  In 
the  case  of  the  400  Watt  metal  halide  lamp  on  the  billboard,  about  half  of  the  en¬ 
ergy  used  by  the  lamp  is  disipated  by  heat  [5].  The  spectral  distribution  functions 
used  by  the  model  are  shown  in  figure  4. 
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Figure  4:  UV  and  IR  normalized  spectral  output  of  Metal  Halide  lamp 

and  enclosure.  The  IR  spectra  is  a  gray  body  spectra  computed  from  bulb  and 
lamp  temperature.  The  UV  spectra  is  from  a  Sylvania  Metalarc^^  digitized  from 
Sliney  [4] 


4  Aircraft  Signatures 

For  aircraft  mounted  UV  sensors,  afterburners  are  of  special  interest  for 
two  reasons.  First,  they  can  generate  a  large  quantity  of  UV  radiation  which  can 
interfere  with  the  operation  on  board  sensors.  Second,  as  fast  moving  combustion 
sources,  they  can  cause  false  positives. 

Modeling  in  this  area  started  with  extensive  ground  measurements  of  F- 
16  and  F-15  aircraft  at  Eglin  AFB  [7]  [6].  These  measurements  included  corre- 


lated  IR  and  UV  measurements  but  only  consider  the  aircraft  operating  tied  to  the 
ground  at  sea  level.  Additional  fixed  wing  measurements  at  lOKft  were  made  of 
afterburners  and  countermeasure  flares  using  a  CMWS  prototype  [8].  While  these 
measurements  were  not  well  calibrated,  they  include  TSPI  and  provide  some  use¬ 
ful  bounds  on  UV  signatures. 

Because  the  data  sets  used  for  the  PEAS  afterburner  models  are  limited, 
it  is  useful  to  compare  the  output  of  radiative  plume  models  like  SPF-III[9]  / 
SPURC[10].  If  the  model  predictions  agree  with  the  measured  signatures,  they 
can  be  used  to  provide  guidance  for  the  scaling  of  plume  signatures  to  the  opera¬ 
tional  velocities  and  altitudes  not  measured  during  the  ground  test.  The  compar¬ 
ison  is  also  useful  because  the  infrared  signatures  in  the  simulaiton  were  already 
dependent  on  radiative  plume  models.  If  the  model  could  use  the  same  inputs  to 
accurately  predict  the  plume  signature  in  both  regions,  confidence  in  the  model 
would  be  increased. 

Since  SPF-III  and  SPURC  had  not  been  utilized  to  model  aircraft  after¬ 
burner  plumes  at  the  wavelengths  of  interest,  the  capabilities  and  limitations  of 
the  models  were  evaluated  to  provide  information  on  the  sensitivity  of  the  models 
to  changes  in  parameters  and  inputs. 

Our  approach  was  to  compare  the  predicted  spectra  and  plume  size  us¬ 
ing  the  same  inputs  for  both  the  infrared  and  UV  portions  of  the  spectrum  to  the 
ground  data  collected  at  Eglin.  The  details  of  the  code  analysis  are  beyond  the 
scope  of  this  paper,  but  their  was  good  agreement  of  model  predictions  to  mea¬ 
sured  infrared  and  UV  plume  lengths  and  overall  radiant  intensities.  This  agree¬ 
ment  required  some  changes  to  the  turbulence  model  in  SPF-III  and  some  changes 
to  the  plume  chemestry  model  in  SPURC  [11].  The  infrared  spectra  also  agreed 
well,  but  the  UV  spectra  did  not  as  the  model  predicts  a  color  temperature  that 
was  significantly  lower  than  the  data  suggest. 

The  primary  source  of  UV  photons  in  afterburning  plumes  is  chemilumi¬ 
nescent  emission  from  CO  -i-  O  resulting  from  the  oxidation  of  soot  in  the  plume. 
The  code  was  traced  to  fundamental  CO  -i-  O  emission  properties  reported  in  [12] 
and  the  code  was  found  to  be  making  predictions  consistant  with  this  data. 

Based  on  the  good  match  to  UV  radiant  intensity  and  plume  size,  the  UV 
output  from  the  model  was  used  formulate  a  tables  for  altitude  and  velocity  scale- 
ing  functions  P{a)  and  y(u).  The  spectral  distribution  and  aspect  angle  factors 
were  retained  from  measured  data.  Afterburner  signatures  with  steady  state  and 
start  up  transients  are  included  in  the  database. 


5  Conclusion 

This  paper  outines  the  development  of  a  UV  signature  database  of  emitters 
that  match  IR  emitters  in  a  closed  loop  simulation. 
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